home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / src / exampleCode / MP / kernel_switch / mpy_mips2.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.5 KB  |  68 lines
  1. c 
  2. c Copyright (C) 1993, Silicon Graphics, Inc.
  3. c All Rights Reserved.
  4. c
  5.     subroutine mpy_mips2(a,lda,b,ldb,c,ldc,m,n,p,alpha)
  6. c
  7.     implicit double precision(a-h,o-z)
  8.     real*8 alpha,beta
  9.     dimension a(lda,p),b(ldb,n),c(ldc,n)
  10.     integer p
  11. c
  12.     double precision t
  13. c
  14.     do  100 j=1,(n-1),2
  15.         do 50 i=1,(m-3),4
  16.             t00 = c(i,j)
  17.             t10 = c(i+1,j)
  18.             t20 = c(i+2,j)
  19.             t30 = c(i+3,j)
  20.             t01 = c(i,j+1)
  21.             t11 = c(i+1,j+1)
  22.             t21 = c(i+2,j+1)
  23.             t31 = c(i+3,j+1)
  24.             do 30 k=1,p
  25.             t00    = t00    + a(i,k)*b(k,j)
  26.             t01    = t01    + a(i,k)*b(k,j+1)
  27.             t10    = t10    + a(i+1,k)*b(k,j)
  28.             t11    = t11    + a(i+1,k)*b(k,j+1)
  29.             t20    = t20    + a(i+2,k)*b(k,j)
  30.             t21    = t21    + a(i+2,k)*b(k,j+1)
  31.             t30    = t30    + a(i+3,k)*b(k,j)
  32.             t31    = t31    + a(i+3,k)*b(k,j+1)
  33.  30    continue
  34.                         c(i,j) = t00
  35.                         c(i,j+1) = t01
  36.                         c(i+1,j) = t10
  37.                         c(i+1,j+1) = t11
  38.                         c(i+2,j) = t20
  39.                         c(i+2,j+1) = t21
  40.                         c(i+3,j) = t30
  41.                         c(i+3,j+1) = t31
  42.  50    continue
  43.         do 80 i=i,m
  44.         t00 = c(i,j)
  45.         t01 = c(i,j+1)
  46.         do k =1,p
  47.             t00    = t00    + a(i,k)*b(k,j)
  48.             t01    = t01    + a(i,k)*b(k,j+1)
  49.         enddo
  50.         c(i,j) = t00
  51.         c(i,j) = t01
  52.  80     continue
  53.     
  54.  100     continue
  55.  
  56.         do j=j,n
  57.                 do k=1,p
  58.                         t00 = b(k,j)
  59.                         do i=1,m
  60.                                 c(i,j) = c(i,j) + a(i,k)*t00
  61.                         enddo
  62.                 enddo
  63.         enddo
  64.  
  65.  
  66.     return
  67.     end
  68.